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In recent years, there has been renewed interest in line-start synchronous 
reluctance machines (LSSRMs) due to their simple construction, 
magnet-free rotor, and low cost. To improve control performance, design 
optimization, and fault diagnostic analysis of these machines, it requires 
accurate estimation of their electromagnetic characteristics using detailed 
and time-consuming finite element analyses (FEAs). In this paper, 
inductances and electromagnetic torque of the LSSRM were calculated using 
the combination of winding function analysis and conformal mapping 
instead of FEA. The hybrid approach can be applied to the prediction of 
motor behavior, taking into account all space harmonics of the air-gap 
permeance without any restriction on the rotor saliencies and topologies. The 
influence of the core saturation, stator slots, and rotor bars were also 
considered. The results obtained by simulations were compared with FEA in 
terms of accuracy and computational time. 
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1. INTRODUCTION 

Over the past decade, there has been sustained research activity on the line start synchronous 
reluctance machines [1]. Several studies have shown that these types of machines and their online starting 
ability are suitable for a wide variety of applications and are particularly attractive for variable-speed 
applications such as industrial conveyors [2], [3]. The line-start synchronous reluctance machine (LSSRM) 
starts as an asynchronous machine based on the induction-motor principle. Then, the machine is able to 
operate as a synchronous motor, reaching synchronous speed thanks to the reluctance torque [4]. In order to 
analyze this type of machine, several analytical and numerical approaches have been developed [5]. 
However, due to the special topology of the rotor, which combines saliency and rotor cage slots under 
magnetic saturation conditions, most of the proposed analytical approaches are rarely able to offer fast and 
accurate analyses of self-and mutual inductance characteristics [6]. The LSSRM electromagnetic torque splits 
into two components: a reluctance pulsating component and a constant component which is obtained when 
the machine is fed by a sinusoidal current [6]. Electromagnetic torque and machine armature winding 
inductances are both functions of rotor position and phase current. The torque and winding inductances can 
both be well determined by two well-known methods, numerically with finite element analysis or analytically 
with the classical winding function theory [6], [7]. An accurate result can still be obtained using finite 
element analysis (FEA). However, this approach is often difficult to apply and leads to significant 
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computational time [8]. The second method performs the integral over the mechanical position of the ratio of 
stator winding turns distribution to the air-gap permeance waveform [9], [10]. This method is highly effective 
because of its simplicity, but has insufficient accuracy due to the air-gap permeance of Fourier’s function. 
This function, which is the sum of superimposed sine and cosine waves of changing amplitude, frequency, 
and phase, represents the decomposition of the air-gap permeance variation caused by the rotor saliency 
combined with the stator and rotor slot effects. The saturation influence was estimated using an equivalent 
air-gap length. Under transient conditions, each part of the motor core experiences different magnetic flux 
densities and different flux path lengths. Based on the equivalent air-gap length technique, the substituted 
air-gap length was computed and adjusted to the saturated conditions of the machine, which implied its 
enlargement. The LSSRM was described almost completely in terms of the geometrical characteristics of its 
air-gap, such as the air-gap length, slot opening width, and the number of slots. Once the air-gap geometry 
was defined, the execution of the proposed hybrid winding function theory (HWFT) developed by [11] was 
used to complete the self and mutual inductance calculation process. At several points, the approach does not 
require any dealing with complex physics principles involved in the whole LSSRM device structure. The 
technique closely follows the performance of the traditional FEA in terms of accuracy and computational 
time. In this work, two methods were used and compared for inductances and electromagnetic torque 
calculation, FEA and the new HWFT method using conformal mapping for the complex transformation of 
the air-gap permeance function. The obtained rotor permeance function was directly used in the analytical 
expression to compute self and mutual inductances. Both HWFT and FEA were compared in terms of 
accuracy and time-consuming process. 


2. MATERIAL AND METHOD 
2.1. Description of the LSSRM model 

The rotor geometry of the LSSRM is similar to the rotor of an induction motor with a squirrel cage. 
The shape of the rotor is modified to a salient rotor by cropping four parts from the cylinder periphery core, 
which represents the four-rotor poles. The cropped areas around the rotor periphery were limited to 
one-fourth of the rotor radius to ensure a low ratio between d-axis and q-axis rotor inductance. Figure 1 
shows the cross-section of the LSSRM with the standard three-AC stator structure with 40 slots. Table 1 
shows the detailed machine dimensions. The study takes into account the saturation effect. 
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Figure 1. Line-start synchronous reluctance motor cross section 


Table 1. Dimensions of the machine 


Quantity and Symbol Value 
Stator outer diameter (D) 160 mm 
Inner radius of the stator (R) 45.5 mm 
Outer radius of the rotor (r) 45 mm 
Constant air-gap length (e) 0.5 mm 
Active axial length (L) 70 mm 
Pole arc/pole pitch (£) 0.5 
Number of poles pairs (P) 2 (4 poles) 
Number of slots in the stator (n) 36 
Turns per slot (N) 115 
Stator slot pitch (T) 10 degrees 
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2.2. Hybrid winding function theory and conformal mapping 

The classical winding function theory was of great help for the analysis of standard electric 
machines, such as induction motors and DC motors. This method is still useful for transient and steady-state 
analyses of salient machine structures such as salient pole generators, switched reluctance motors, and 
permanent magnet synchronous generators [12]-[16]. The winding function expression consists of a definite 
integral involving the geometrical position of the phase winding which represents the magnetomotive force 
distribution across the air-gap obtained by (1) and the air-gap permeance profile that includes slotting and 
saliency effects [17], [18]. The geometrical position of the phase winding N, can be expressed by the Fourier 
expansion as (1): 


Ny) = DR, A(2h — 1) cos (Qh -1)(¢- =) i 


2 


where A represents the numerical Fourier expansion of the total winding factor, which is equal here to the 
product of the distribution factor, the skewing, and the pitch factor [19]. The symbol k represents the stator 
phase a, b, or c. The waveform of N, obtained by (1) for the phase a, b and c respectively is shown in Figure 2. 

However, the classical winding function theory analysis gives only accurate results for machines 
with small smoothed air gap topologies. The rotor permeance equivalent function using the Fourier series 
often does not provide an accurate model for slotted structures [20], [21]. 
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Figure 2. Winding function of the stator phases a, b and c 


2.3. Rotor permeance function using conformal mapping 

The technique proposed by [11]-[13], reviewed the air-gap permeance function using the MATLAB 
Schwartz-Christoffel (SC) Toolbox for conformal mapping transformation [14], [22]. The authors propose a 
numerical identification of the conformal map using the complex transformation technique instead of an 
explicit Fourier expansion for an accurate profile of the permeance. This was performed using a logarithmic 
complex function depicted in (2), which converts the polar coordinate of the circular real air-gap geometry 
named the S-plane into a complex plane in the Z-plane. 


Z = log(S) (2) 


Figure 3 shows the air-gap geometry in Figure 3(a) and its permeance identification after SC mapping 
transformation in Figure 3(b). The numerical SC Toolbox [23] was used to obtain an open complex 
configuration of the planar polygon geometry which represent the map of the combined stator slots and rotor 
saliency with cage slots. 


Z= As” m w= w,)™ tdw + cst 8) 


Where A and cst are the equation constants relative to the integration, n is the number of vertices in the 
z-plane, W,,...,W,, are the points located in the canonical rectangle in the w-plane corresponding to the 
vertices and a, are the angles of the polygon. 
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Figure 3. The inverse air-gap permeance mapping (a) physical air-gap and (b) flow chart of the inverse air- 
gap permeance identification of LSSRM by the usage of SC mapping 


Using expression 3, the transformation leads to a canonical smooth rectangle called the w-plane of 
the stator geometry configuration, keeping the air-gap permeance easy to determine. The canonical rectangle 
represented in the w-plane provides the corresponding inverse gap waveform as shown in Figure 4. 
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Rotor position (Elec .deg.) 
Figure 4. Saturated and unsaturated inverse air-gap waveform 

The curve shows the permeance waveform for the unsaturated core material in a dashed line versus 
the saturated core material's. Once the permeance function P,,. is identified by the SC toolbox, the numerical 
Fourier expansion of the inverse gap waveform is calculated and presented by (4). 

Pop (x) = Dna P,cos(h x Xr) (4) 

In order to calculate the flux density distribution along the air-gap circumference, which is produced 
by the stator AC rotating field, the winding function N,(@)is multiplied by the inverse air-gap P, function 
and the phase current flowing in that winding [24]. 


Ba (xr) = UoPsr(@ — Xr) Na(@)ia (5) 


Where ig is the RMS current flowing in stator phase (A), Na The number of coil turns. The magnetic flux 
between two stator phases k and k’ can be expressed as (6). 


Arr (xr) = ia * EE? x) (6) 


Lkr (xr) is the expression of two stator phases k, k' mutual inductance written as (7): 


Lyr (xr) = EE 4 129 (x) (7) 
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where 1480 which is independent from the rotor position [25] represents the term accounting for magnetic 


flux leakage with analytical validation in [26] and the term LED is the magnetizing flux between two stator 


phases. The term a can be expressed using: 


Lpr (Xp) = RLMg Jo” PoelO — Xp) NC) * Neb) dp (8) 


Fourier expansion of (1) and (5) can be substituted into 6 giving: 


Le (X-) = RLU KER- A(2m — 1)cos (com -1) ( - =n(k - »)) Yr, A(2n — 1) cos (cen -—1) ( - 
Enae- 1))) + Bia acos( WC =x) dd (9) 


Which can be written as: 


Lyae(%p) = >) 5 2 Ph* A(2m —1) 


nis 2(k — 1) 
D) 


== )de (10) 


* A(2n — D|” cos(h(g — xr)) * cos( (2m — 1) (p — 


< 1). 


* cos( (2n — 1)(ọ — 


where, [lg is the air-gap permeability, i, the current of the phase (A), x, and @, are the angular position of the 
rotor respect to the location of the stator reference axes and the position along the inner stator region 
respectively. The rotor outer radius of the LSSRM is a function of the angular position x,, and the inverse 
air-gap length, presented by the expression 2 is also a function of the rotor position x, [27]. 

In other words, the salient air-gap permeance varies and alternates periodically as the rotor rotates. 
As a result, two specific paths of the stator magnetic flux are identified. The first one is aligned with the 
direct rotor axis representing the pole's face region, which is standardly related to the small air-gap, and the 
second one is aligned with the rotor quadrature-axis position, which is related to the large air gap. Figure 5 
shows the two-flux density; one is a high permeability path, the flux lines follow rotor iron paths crossing the 
minimum air-gap length, the second represents the maximum air-gap length with a low permeability path. 

The direct and quadrature calculation of the magnetic flux density was done using one of these 
following calculation and curve arrangement: Shifting by an angle of 45° the spatial angle of the rotor 
permeance curve with the respect to the magnetomotive force in an electric degree. Or, shifting with the same 
angle the magnetomotive force of the three stator phases with the respect of the rotor permeance. 


2.4. Stator tooth saturation modeling 

The approach described in [11] was found to be an ideal framework for the calculation of the 
magnetic saturation of various parts of the core. The method can be easily adapted to tooth saturation 
analysis. Besides dividing the stator tooth only into a maximum of two regions, in our approach, the stator 
tooth body is represented by several polygons. Each polygon represents a region with different lengths and 
areas. The stator tooth is shown in Figure 6, where regions are delimited by different length polygon 
boundaries on the surface with crosshatch. By definition, the winding function theory (WFT) assumes the 
relative permeability of the LSSRM core is infinite. In this study, the saturation effect is taken into 
consideration using the technique proposed by [11]. 

The authors use the WFT method by means of infinite permeability assumptions, in which the 
saturation is substituted by adjusting the equivalent air-gap length. A suitable value for the length of the air 
gap has thus to be determined for all stator teeth. The saturation effect is simulated by increasing the air gap 
length between the stator teeth and the rotor periphery. The longer the air-gap, the more saturated the core. 
Saturation modeling includes the following steps: i) the virtual coil is wound so that it has a one-turn coil on 
each tooth; ii) the incremental self and mutual inductance were carried out depending on the rotor position; 
and iii) the current through each virtual coil must be determined. The tooth’s MMF is equivalent to the 
current of it is corresponding tooth. At each rotor position the flux linkage for the corresponding virtual tooth 
coil are calculated. 

The equivalent length deltaLg is the sum over the number of tooth sections ranged from 1 to the 
maximum section number Nmax of the length body over the body sections times their respective 
permeability. here we choose only one section length due to the fact that the core permeability is not 
changing. The equivalent air-gap length deltaLg is presented by the following generalized (11): 
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At poc nmax Ityoay) 
deltaLg = —#2¢ + lro” Du baay 


Hr, n=1 At poay() 7 Lt vody (n) 


shoe 


(1) 


where the parameters Itpoay and Itshoe represents the body and the shoe length, Atpoay and Atshoe represents the 
cross area for the body and the shoe respectively and uroay and prtshoe for the relative permeability. Although 
the Schwartz-Christoffel Toolbox is a useful tool, it suffers from performance limitations due to the rotor slot 
opening. 

The Schwartz-Christoffel MATLAB toolbox may not always be able to accurately calculate the 
equivalent polygon which contains the entire stator and rotor slotted geometry. For the studied model, the 
stator and rotor slotted geometry design was limited to only the slot opening region, which represents a third 
of the total slot height. The diagnostic algorithm process may lead to inaccuracies because of these 
incomplete slot geometries, which are not always considered effectively. 
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Figure 5. Air-gap flux density in d-axis and q-axis Figure 6. Stator tooth with virtual coil 


3. RESULTS 

The self and mutual inductances are calculated using expression 8 for various rotor positions; 
Figure 7 clearly shows ripples on the inductance profiles due to the slotting effects. the figure of the 
inductance waveforms shows the rotor saliency influence for various rotor positions that leads to an 
additional pulsating torque and supplements the average torque, when the rotor saliency ratio is high the 
torque profile will distort and may lead to a large a pulsating torque component. 


3.1. Analytical calculation of the air-gap electromagnetic torque 
The air-gap electromagnetic torque TX, of the kt” armature phase, is calculated using the following 
derivative of the magnetic co-energy expression: 


(k) 
TemW = [eee (12) 
i=Cst 


Oxy, 


The total torque is represented by the sum of the three armature phases contribution. 


Lis 
g OLC,) y i 


Fz = 1 ee i i l. (13) 
airgap n. em 2 [ as bs cl a 
les 


Where Tx is the electromagnetic torque contribution of each armature phases and the total energy stored 
W.o-energy tepresents the co-energy which is valid only for linear systems. ias, ipẹ and i., are the three-phase 


current supplying the stator winding. œ is the angle that must be regulated to maximize the electromagnetic 
torque. here the control angle œ was set to a value of 0.3245 rad. 
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igs = y 2İÍm sin (xr = =) 
$ 7 j T 
ibs = y 2im sin (x, =o ((4+ a)r)) (14) 
igs = 4 Zim sin (x, = : ((2 + a)n)) 
Lsaa Lsab Lsac 
L(x,) = | Lsba Lspp  Lsbe (15) 
Lsca Lscp Lsce 


Where L(x,.) represents the inductance matrix. this matrix allows the air-gap electromagnetic torque 
computation by means of (16). The current vector is shifted from the phase reference “a” with 45° electrical). 
The torque characteristics versus rotor position in electrical angle can be calculated as it is shown in Figure 8. 
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Figure 7. Winding function analysis Figure 8. Calculated torque versus rotor position 
with HWFT (5=45°, Irms=1 A); hybrid winding function method 


3.2. Finite element analysis and hybrid winding function analysis comparison 
3.2.1. The flux density distribution in the air-gap 

Two-dimensional (2D) magnetostatic finite element analysis was carried out with the same set of 
parameters used for the winding function analysis. The 2D FEA was performed on the basis of a commercial 
software package [28]. Figure 9 (a)-(d) shows the radial air-gap flux density with respect to electrical degree 
on the d-axis in Figures 9(a) and 9(c) and q-axis in Figures 9(b) and 9(d) respectively. The flux density of the 
direct axis is slightly lower than the flux density of the quadrature axis, since the flux distribution in the d 
axis is denser than the quadratic axis flux distribution. The FEA results revealed good agreement with the 
HWFT results. 


3.3. The finite element computation of stator inductances 

In the stator stack, the phase winding self-inductance "a" can be calculated by using numerical 
integration techniques as well as the estimated value of the co-energy stored in the magnetic field that 
couples the adjacent phases. The inductance and the co-energy are expressed by (14) and (15): 


W, 
Lie = (16) 
Wert = $ Bp Hy dQ (17) 


where, W,,,’ is the magnetic field energy located in the air-gap linking conductor k with conductor k’, I, 
represents the current in the "k" conductor, B,is the magnetic flux density where the current flowing in the 
conductor k is equal to one ampere, Hp, is the magnetic field where the current flowing in the conductor k’ is 
also equal to one ampere. For multi-turn conductors, the inductance quantity is given by (16): 


Lnet = N?Lmatrix (18) 
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where Lmatrix represents the matrix with self and mutual components for the three stator phases. The results 
obtained for self and mutual inductance profile by FEA (compared to those obtained with the HWFT) are 
shown in Figure 10. The self and mutual inductance harmonics are depicted respectively in Figure 11(a) and 
Figure 11(b). Both analysis of the finite element method and winding function theory expose similar results 
except for the direct component of the self-inductance [29]. The difference in the direct component is due to 
the leakage flux in the stator slots. The leakage flux is neglected while developing winding function 
calculations. 

Computational time required for finite element analysis at a numerical resolution of 2/100 rad takes 
around 11 minutes for a Xenon processor work station computer running Windows 7 with 12 GB RAM). The 
computational costs are in the range of 1 minute using the same PC with the winding function theory coded 
in the MATLAB programming language. Figure 12 shows the spectrum for the air-gap electromagnetic 
torque with respect to the amplitude over a frequency range of 0 to 40*f Hertz. 
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Figure 9. Comparing equipotential of (a), (b) the air-gap flux density with respect to electrical degree on the 
d-axis and (c), (d) equipotential in q-axis, respectively 
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Figure 10. Winding function analysis vs FEM: self-inductance curve profile of the phase "a"; mutual 
inductance curve profile between a and b 
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3.4. The electromagnetic air-gap torque computation 
The numerical approach uses the virtual work principle to compute the torque on the rotor region. 


Across the air-gap, the electromagnetic torque expression is given by the following relationship expressed by 
(17): 


Tia 5 l F Bp) dH (p)dv (19) 


where B represents the flux density and H is the air-gap flux intensity. Both B and H are a function of the 
angle p. When computing electromagnetic torque, the mesh was mostly refined within the air-gap region, 
however, in order to save a lot of the computation time, the mesh operation performed with the dedicated 
software is adjusted to keep the mesh node number as low as possible during the iterations [30]. As expected 
from the inductance waveforms results, the torque profile shows a ripple caused by the stator and rotor 
slotting effect. The electromagnetic torque result illustrated in Figure 13 is in accordance with the result 
obtained with HWFT. The computational time required for finite element analysis takes around 660 seconds 
at a resolution of 7/100 rad, the same processor computes the inductance profiles with the same resolution 
within 90 seconds. 
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Figure 11. Harmonic terms: (a) self-inductance, (b) mutual inductance for the two methods 
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(5 = 45°, Ipms=1 A); both HWFT and FEA 


4. CONCLUSION 

The main objective of this work is to investigate the hybrid winding function theory based on WFT 
and conformal mapping for improving performance and analysis of line start synchronous reluctance motors. 
The combination of the two techniques allows for accurate calculation of all space harmonics associated with 
the air-gap permeance without any restriction on the rotor saliencies, stator slots, and core saturation. The 
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results of calculations using hybrid winding function theory agree quite reasonably with the FEA results. This 
approach is well-suited for fast geometric parameter treatment when the LSSRM models are analyzed under 
saturated conditions. The hybrid winding function theory can also be used for design, optimization, and fault 
diagnosis in electric machines. 
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